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Abstract 

In this paper we use convex analysis and variational inequality methods 
to establish an existence result for a model of infinitesimal rate-independent 
gradient plasticity with kinematic hardening and plastic spin, in which the lo¬ 
cal backstress tensor remains symmetric. The model features a defect energy 
contribution which is quadratic in the dislocation density tensor Curl p, giving 
rise to nonlocal non-symmetric kinematic hardening. Use is made of a recently 
established Korn’s type inequality for incompatible tensor fields. The solution 
space for the non-symmetric plastic distortion is naturally H(Curl) together 
with suitable tangential boundary conditions on the plastic distortion. Con¬ 
nections to other models are established as well. 
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1 Introduction 

In the past twenty years, there have been several experimental investigations for metallic 
and ceramic materials which show that the elastic-plastic deformation of those materials 
are size-dependent for sufficiently small scales. This phenomenon cannot be predicted 
by conventional theories of plasticity, which do not include any material length scales. 
Hence, there clearly appeared a gap between micro-mechanical plasticity and classical 
continuum plasticity. The purpose of the enhanced gradient plasticity theories is to 
formulate a constitutive framework on the continuum level which is used to bridge 
that gap. In this paper we will only discuss phenomenological models of isotropic 
polycrystalline plasticity excluding the important case of single crystal plasticity. 

There is now an abundant literature, and research activities towards the develop¬ 
ment of models of gradient plasticity which capture better the observed size-dependency 
mentioned above. In the works of Mfihlhaus and Aifantis [33] and Gudmundson [20], 
the yield-stress is set to depend also on some derivative of a scalar measure of the accu¬ 
mulated plastic distortion. In the works of Gurtin and Anand [21], Gudmundson [20] 
and Neff et al. [34], the yield-stress is not modified, but the free-energy is augmented 
by a term involving the dislocation density. Also, it is assumed in [21] that the plastic 
flow is governed not necessarily by the stress deviator (as in classical plasticity), but 
more generally by microstress tensors that also satisfy a balance law. 

In [33, 20, 21], the plastic distortion variable is assumed to be symmetric. However, 
as Gurtin and Anand note [21, p. 1626]: ’’...unless the plastic spin is (explicitly) 
constrained to zero, constitutive dependencies on the Burger tensor necessarily involve 
dependencies on the (infinitesimal) plastic rotation”. Note that even in classical plastic¬ 
ity, the effect of plastic spin has been studied by several authors like Dafalias [11, 12], 
Mandel [28, 29] and Kratochvil [25, 26], who were the first to suggest that a com¬ 
plete macroscopic plasticity theory must include constitutive relations involving also 
the plastic spin. 

On the other hand, though there are several theories of gradient plasticity available 
in the literature, the results of mathematical analysis for these problems are still rather 
scarce. The first result of mathematical analysis on a model of gradient plasticity 
was due to Djoko et al. [14]. While the developments by Gurtin-Anand [21] were 
done for viscoplastic bodies, the well-posedness of that model is considered by Reddy 
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at al. [46] for the rate-independent problem with isotropic hardening and with both 
energetic and dissipative length scales involved. Also, computational aspects of the 
model, based on the work [23], are studied in [7] and are devoted exclusively to single 
crystal plasticity. Let us mention that, the purely energetic version of the Gurtin-Anand 
model, i.e., when the dissipative length scale £ = 0, is not yet treated but will be done 
in this paper via the identification with the irrotational version of Ebobisse-Neff [17] 
presented in Paragraph 3.4. We also study in Section 5, the purely energetic case of the 
Gurtin-Anand model with linear kinematic hardening. Another existence result for the 
rate-independent problem of the Gurtin-Anand model was obtained by Giaccomini and 
Lussardi [19] within the energetic-approach developed by Mielke [31, 32] and it has also 
been proved that the model converges in a suitable sense to a formulation of classical 
perfect plasticity proposed in [13] whenever the energetic and dissipative length scales 
go to zero. 

Neff et al. proposed in [34] a model of finite strain gradient plasticity based on 
the multiplicative decomposition including phenomenological Prager type symmetric 
linear kinematical hardening and nonlocal kinematical hardening due to dislocations. 
The model is from the outset non spin-free (the plastic distortion p is not symmetric) 
and its linearization leads to a thermodynamically admissible model of infinitesimal 
plasticity involving only the Gurl of the non-symmetric plastic distortion p. The well- 
posedness of the linear model is addressed as well, when formulated as a variational 
inequality. 

In [17], we have studied the well-posedness of this model described within the frame¬ 
work of the dual formulation for isotropic hardening by 


Additive split of distortion: 

Vu = e + p 

Additive split of strain: 

sym Vtt = e = e® + , e® = sym e 

Equilibrium: 

Div cj -|- / = 0 with a = C.e® 

Dissipation inequality: 

J [(fj — ^ Gurl Gurl p, p) — k 2 'y jjdx > 0 

p € dX{a — p L'lCmlCurlp), f = \p\ 

X is the indicator function of the set 
of admissible stresses 

Flow law in dual form: 


Table 1: The model with plastic spin and isotropic hardening in [17]. 


In this paper, we settle some of the open questions raised in [17]. In that paper, 
we dealt with the isotropic hardening case only and here we would like to extend our 
analysis to the local linear Prager type kinematical hardening model which, in the dual 
formulation of classical plasticity has the flow law 

€ dX(a — b) and b = pki£^, (1.1) 

in which b is the symmetric backstress tensor. Notice that ( 1 . 1)2 can be explicitly 
integrated (with proper initial conditions e^(0) = 0, 6(0) = 0) to yield b = pkie^ and. 
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substituted in (l.l)i, gives 

G dX(cr — fiki e^). 

Applying the same reasoning to our nonlocal gradient model, we want to add a local 
backstress similar to ( 1 . 1 ) 2 . Hence, we define 

b = fiki symp , 

in which we conserve the symmetry of the local backstress tensor b. So, we get similarly, 
p G dX{a — b — CurlCurlp) and b = pki symp 

and integrating as above yields 
d V 

p G clX(S£;) C with = O'— p/ci symp — Curl Curlp. (1.2) 

I dev 2 j e \ 


Here, is the elastic Eshelby tensor driving the plastic evolution. Note immediately 
that it is only the nonlocal backstress contribution p Curl Curl p which is responsible 
for the appearance of plastic spin or not. In order to substantiate this claim, set Lc = 0 
in ( 1 . 2)2 and consider 


dev Ee 
^ |dev Ee\ 


G Sym(3). 


Assuming p(0) G Sym(3), we get that p{t) G Sym(3) for every t and hence we may 
replace p with := sym p. 

From a purely mathematical point of view, we could also consider the case of a 
non-symmetric local backstress tensor, i.e.. 


p G 9X(cJ — 6 — CurlCurlp), and b = pkip, 
which integrates to 

p G fix (cr — p A:i p — p Curl Curl p). (1.3) 

In this case, the following mathematical analysis would be drastically simplified since 
no new estimate of the Korn type on incompatible tensor fields (see [40, 37, 38, 39]) is 
needed. The solution space is then trivially H(Curl). Notice that (1.3) also still reduces 
to a formulation of classical plasticity in terms of a symmetric plastic strain tensor 
if the energetic length scale Lc vanishes and the initial plastic distortion p(0) is chosen 
to be symmetric. To see this, consider for Lc = 0, the equation 

^ dev a — pkip 
^ I dev a — pkip\ 


The format of the equation, as far as classical solutions is concerned, is of the type 

p = S{t) - ap{t), p(0) = po G Sym(3), (1.4) 
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where S(t) G Sym(3) and a G M can be assumed given. Clearly, the system (1.4) has 
only symmetric solutions p{t). 

The total energy would be of the type (isotropic elastic response for simplicity) 
/i|sym(Vii-p)p + ^ |tr(Vu-p)p + ^ |Curlpp , 

immediate //(Curl )-control of p 

which is, however, not invariant w.r.t. the transformations 

u ^ u + Ax + b, p ^ p A , 

for constant skew-symmetric A G so(3) and constant translation 6 G which rep¬ 
resent superposed Euclidean motions on both the displacement and plastic distortion. 
Therefore, the choice symp in the backstress evolution is mandatory by Euclidean in¬ 
variance; the linear kinematic hardening must be based on a symmetric backstress 
tensor. 

Before we present our analysis of the model with linear kinematical hardening and 
plastic spin, we find it important to first present, using the convex analytical setting, 
a summary of those few models of infinitesimal gradient plasticity in the literature for 
which a mathematical analysis is now available. Precisely, we present in Section 3 the 
model by Miihlhaus-Aifantis [33] as analyzed in [14, 15], the model by Gurtin-Anand 
[21] as studied in [46], the models with plastic spin analyzed in [34, 17] and their 
irrotational version in [41] and we highlight some interconnnections. 

Let us first fix some notations and definitions which will make the paper more clear 
and readable. 

2 Some notational agreements and definitions 

Let n be a bounded domain in with Lipschitz continuous boundary <911, which is 
occupied by an elastoplastic body in its undeformed configuration. Let T be a smooth 
subset of 50 with non-vanishing 2-dimensional Hausdorff measure. A material point 
in 0 is denoted by x and the time domain under consideration is the interval [0,T]. 
For every a, 6 G M^, we let {a,b)^3 denote the scalar product on with associated 
vector norm |a |^3 = («)tt)][^3. We denote by R^^^ the set of real 3x3 tensors. The 
standard Euclidean scalar product on R^^^ is given by {A, B)-^x3 = tr , where 

denotes the transpose tensor of B. Thus, the Frobenius tensor norm is |A|^ = 
{A, A)j^3x 3. In the following we omit the subscripts R^ and R^^^. The identity tensor 
on R^^^ will be denoted by 1, so that tr(A) = (A, 1). We let Sym(3) ;= {X G 
j^3x3 I (Jenote the set of symmetric tensors, the Lie-Algebras so(3) ;= {X G 

]^3x3 I _ —X} of skew-symmetric tensors and s[(3) ;= {X G R^^^ | tr (X) = 0} 
of traceless tensors. For every X G R^^^, we set sym(X) = ^(^X + skew(X) = 

4(X — X"^) and dev(X) = X — ^tr (X) 1 G s[(3) for the symmetric part, the skew- 
symmetric part and the deviatoric part of X, respectively. 
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The body is assumed to undergo infinitesimal deformations. Its behaviour is gov¬ 
erned by a set of equations and constitutive relations. Below is a list of variables and 
parameters involved in various models of infinitesimal gradient plasticity presented in 
this paper: 

• u the displacement of the macroscopic material points; 

• p the plastic distortion variable is a non-symmetric second order tensor, incapable 
of sustaining volumetric changes; that is, p G sl(3); 

• e = Vrt — p the elastic distortion is a non-symmetric second order tensor; 

• = sym p the symmetric plastic strain tensor; 

• e® = sym {S/u — p) the symmetric elastic strain tensor; 

• a the Cauchy stress tensor is a symmetric second order tensor; 

• cjy the yield stress; 

• / the body force; 

• Curl p = — Curl e the dislocation density tensor; 

• the microstress tensor is a second order deviatoric symmetric tensor; 

• = {niijk) the micro-polar stress tensor is a third order tensor deviatoric sym¬ 
metric in the first two indices i and j. That is, rriijk = rrijik and muk = 0; 

• 7 the accumulated plastic strain. 


For isotropic media, the fourth order elasticity tensor C is given by 

C.X = 2p dev symX -|- k tr(X)l = 2p symX -|- A tr(X)l (2.1) 

for any second-order tensor X, where p and A are the Lame moduli satisfying 

p > 0 and 3A -|- 2|U > 0 , 
and K > 0 is the bulk modulus. 

These conditions suffice for pointwise ellipticity of the elasticity tensor in the sense that 
there exists a constant mo > 0 such that 

(X,C.X) > mo|symX|2. (2.2) 


The space of square integrable functions is L^(I2), while the Sobolev spaces used in this 
paper are: 


H^(I2) = {u € L?‘{VL) I gradrt G , 

ll«ll/fi(o) = ll«lli2(o) + Ilgradw||i2(t^) , 

H(curl; 17) = {u G L‘^{Q) \ curlu G L^(17)} , 

ll^llH(curl;0) = Minn) + I|curlu||i 2 (j^) 


grad = V, 

Vn G H^(17), 
curl = V X , 

, Vu G H(curl; Q). 
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For every X G with rows we use in this paper the definition 

of Curl X in [34, 47]: 

/ curlXi “ “ \ 

CurlX= curlXs-G M^xs ^ 

\ curlXs - - / 

for which Curl Vu = 0 for every v G M^). Notice that the definition of Curl X 

above is such that (Curl X)^a = curl (X^a) for every a G and this clearly corre¬ 
sponds to the transpose of the Curl of a tensor as defined in [21, 22], 

In Paragraph 3.5, we will need an explicit definition of the linear operator L : 
IR3x3x 3 ^ ]^3x3 


So, for 


CurlX = L.VX VX G R3X'^). 

( 0 -03 02 \ 

03 0 —Oi j Gso(3), 

-02 Oi 0 / 

we consider the operator axl : so (3) —)■ through 

axl(y4) ;= (oi, 02 , 03 )^ , X. u = (axlX) x u, VuGIE 

1 3 _ 1 3 _ 

(axlyl)^ — ^ ^ ^ ^ , 

i,j=l i,j=l 

where Cijk is the totally antisymmetric third order permutation tensor. 
Hence, for every A G 


(2.3) 


1 3 1 3 43 

(axl skew yl)/j — ^ ^ e^jj skew(yl)jj — — ^ ^ c^ijAji — ^ ^ CkijAji 

i,j=^ 

1 ^ 

“ 2 ^ 


*d=i 




■ 


i,j=^ 


Recalling that (curl v)k = ^kijVj,i for every v G C'3(H,r 3[)^ it follows that 

id=i 


1 .3 ^ \ 

(axlskew Vu)fc = ^ ^ = -(curl v)k ■ 

Therefore, we may rewrite 


*d=i 


2 axl skew VXi_ 

Curl X = I 2 axl skew VX 2 ___ ] = L.VX, 
2 axl skew VX 3 _ 
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where L : ]^ 3 x 3 jg gjygj^ }jy 

/ Li.VX ___ \ 

L = (Li,L 2 ,L 3 )^, L.VX = L 2 .VX ___ (2.4) 

\ L 3 .VX ___ ) 

with Lj ; ^ i = 1 , 2 ,3 defined by 

Lj.VX = 2 axlskew VXj. (2-5) 

Hence, we have found the explicit linear operator L ; ^ 3 x 3 x 3 K^^ 3 ^ g^ ^j^^t the 

equality (2.3) holds. 

Notice, that for every X, Y G C^{Q, M^^^) 

3 3 

(Curl X, Curl Y) = ^^(curl Xj,curl Yi) = 4 y~^(axl skew XXj, axlskew Vlj) 

i=\ i=l 

3 3 

= 2 (skew V Xj , skew Vlj) = 2 ^^(skew XXj, VYj). 

i=\ i=\ 

The following function spaces and norm will be used later. 

H(Curl; H, 1^3x3) = {X G L^{n, \ CmlX G L^{n, , 

ll^llH(Curi;n) = m\l2^n) + l|CurlX||2,(^) , vx G H(Curl;H, M^x^), (2.6) 
H(Curl; H, E) := {X : H ^ E | X G H(Curl; H, R^xS)} ^ 

for E := Sym(3), s[(3) or Sym (3) ns[(3). 

We also consider the space 

Ho(Curl; Q, r,R3x3) 

as the completion in the norm in (2.6) of the space {q G C°°{Q, T, R^x3^ | g ^ _ q| _ 

Therefore, this space generalizes the Dirichlet boundary condition 

q X n|r = 0 

to be satisfied by the plastic distortion p or the plastic strain := symp. The space 
Ho(Curl; H, r,E) is defined as in (2.6). 

The divergence operator Div on second order tensor-valued functions is also defined 
row-wise as 

/ divXi \ 

DivX = I divX2 I . 

\ divX3 / 
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3 Some models of infinitesimal gradient plasticity 

3.1 The model by Miihlhaus-Aifantis [33] 

In this model, the yield-stress in the case of isotropic hardening, is set to depend also on 
some derivative of a scalar measure of the accumulated plastic distortion which plays 
the role of the isotropic hardening variable. A summary of the model is presented in 
Table 2. 


Additive split of strain: 
Equilibrium: 

Free energy: 

Yield condition: 
where 

sym Vu = e = e Sym (3) 

Div O' -|- / = 0 with a = C.e® 
l{C.E-,e-) + yk 2 h\^ + yp\Vy\^ 

9 ) = \deva\+ g - ay <0 
g = —/u/c 2 7 -I- 

Dissipation inequality: 

[ [{a,eP) +gf]dx > 0 

Jn 

Dissipation function: 

■D(qr):=l ^yl^l 

' ’ f -|-oo otherwise 

Flow law in primal form: 

Flow law in dual form: 

{a,g) G dV{eP,j). 

7 = A = |e 1 

dev a\ 

KKT conditions: 

Boundary condition on y: 
Function spaces for and y: 

A > 0, 4>{a, g) <0, A (p{a, g) = 0 

7 = 0 on do, 


Table 2 : The model by Miihlhaus-Aifantis [33] as formulated in [14, 15] . 

Under suitable boundary and initial conditions on u, and 7 , the well-posedness as 
well as computational aspects of that model are studied in [14, 15] with the flow law 
formulated in its primal form. 

3.2 The model by Gurtin and Anand [21] as studied in [46] with 
isotropic hardening 

This model is based on the assumption that the power expended by each kinematical 
field be expressible in terms of a system of forces consistent with its own balance. 
Therefore, the model is characterized by two additional stress tensors: a second order 
tensor power conjugate to the symmetric plastic strain and a third order tensor 
power conjugate to the gradient of the plastic strain, which satisfy a microforce balance. 
The latter as well as the equilibrium being derived by the principle of virtual power. 
Since is deviatoric symmetric, it is not restrictive to assume that is deviatoric 
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symmetric and the third order tensor is deviatoric symmetric in the first two indices. 
The model as formulated in [46] is summarized in Table 3 with the purely energetic 
version in Table 6 . 


Additive split of strain: 
Equilibrium: 

Microforce balance: 
where 

sym Vtt = e = e® + G Sym (3) 

Div cj + / = 0 with a = C.e® 

dev a = tP — Div m^, 

tP-. microstress (2"^^ order) 

m^: micropolar stress (3“''^ order) 


Free energy: 

^(C.e®,e®) + ^laLl jCurlePp + ^fik2\'y\" 

) 

Yield condition: 

(j){TP,mP,g) := ^JItpI"^ + + g ■ 

-Uy < 0 

where 

Kiss = mP- mLerg 



(CurleP, CurleP) = (mLerg, Ve^) 


Dissipation inequality: 

[ [{tP, eP) + VeP) + gf] dx > 0, 

Jn 

9 = -hhl 

Dissipation function: 

Viq 0 := [ 

' ’ ' 1^ +00 otherwise 


where 

dP{Q) ■■= \/kP + ^2|Vg|2 


Flow law in primal form: 

irP,ml^„g)€dV{eP,VsP,^) 


Flow law in dual form: 

rP 1 

eP = \ , V£P = \r^ , 1 

~ 9 ~ 9 i 

y = \ = dP{£P) J 

I (•) 

KKT conditions: 

A > 0, (j){TP, 5 -) < 0, A (P{tP, g) = 0 

Boundary conditions fore^: 
Function space for : 

Two length scales: 

= 0 on dD, 

eP{t,-)€H^{n, Sym (3)) 
dissipative i and energetic Lc 



Table 3: The model by Gurtin and Anand [ 21 ] as formulated in [46] with dissipative 
and energetic length scales. 


The well-posedness of the model was studied by Reddy et al. [46]. We would like to 
emphasize here that the starting point of the modelling and analysis is the primal form. 
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The corresponding solution will satisfy the dual form in which it is understood that 
there is an extra consistency condition generated which makes equation (*)2 in Table 
3 possible. The plastic strain variable is assumed from the outset to be symmetric. 
Note that the formulation in [21] as well as in [46] involves the full gradient of the 
plastic strain in the dissipation function, which is controlled in leading then to find 
the plastic strain variable in the Sobolev space Sym(3)) together with the 

possibility to completely prescribe at the boundary. 


3.3 The model with plastic spin in [17] and in [34] 

Unlike the model in [21] with the microstresses and the plastic distortion kept symme¬ 
tric, a model involving the plastic spin is studied in [34] with phenomenological Prager 
type kinematical hardening and in [17] with isotropic hardening. A summary of the 
setting in [17] for the so-called equal spin case is presented in Table 4. 


Additive split of distortion: 
Equilibrium: 

Vu = e+p, e® = syme, = symp 

Div cj -|- / = 0 with a = C.e® 

Free energy: 

- 1 - Curlpp -|- ^pk 2 7 p 

Yield condition: 

(/)(Se, g) := 1 dev T.e\ + g - (Jy < 0 

where 

EE:=a + i = Curl Curl p 


g = - 11^21 

Dissipation inequality: 

[ [{^E,P)+ 97 ] dx>0 

Jn 

Dissipation function: 

1 \q\ if kl < 

’ [ 00 otherwise 

Flow law in primal form: 

{EE,9)edVip,j) 

flow law in dual form: 

. devSe . . , 

*’ = bdevE^|’ T = ^ = l!'l 

KKT conditions: 

A > 0, (l){EE,g)<0, X(l){EE,g) = 0 

Boundary conditions for p: 
Function space for p: 

p X n = 0 on T, (Curl p) x n = 0 on dD \ T 
p{t,-) G H(Curl; Q, 


Table 4; The models with plastic spin in [17] and in [34]. 
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An existence result for the weak formulation of this model is obtained in [17]. The so¬ 
lution space for this model is quite naturally p ^ H (Curl) since the isotropic hardening 
provides an L^-control of the entire plastic distortion p and the energetic defect energy 
adds automatically a control of Curl p £ L‘^. 


3.4 The irrotational version of [17]. 

In [41], the irrotational limit case has been computationally implemented as one of the 
hrst efficient treatments of gradient plasticity. In this model, the plastic distortion p 
remains symmetric and can therefore be written as = symp. A summary of the 
model is presented in Table 5. The well-posedness of this limit case is included in the 
analysis presented in [17]. 

As shown in the next paragraph, this model can also be obtained as a particular 
case of Gurtin-Anand [21] for Z = 0, > 0. Since the dissipative length scale I = 0, 

the solution space is only H (Curl) with the attendant tangential boundary conditions. 
Thus, the existence result in [17] provides also the first existence result for the purely 
energetic Gurtin-Anand model with local isotropic hardening. 


3.5 The Gurtin-Anand model: purely energetic version 

In this section, we would like to compare or find a connection between the model by 
Gurtin-Anand and our irrotational version. To this aim, we consider the defect energy 

1 Curl £P\^ = [L.Ve^’l^ , (3.1) 

where the linear operator L ; —)■ explicitly defined in (2.4)-(2.5) is such 

that 

Curl eP = IL.VeP . 

We recall that 

/ curl £P _\ 

Curl £P = I curl _ ) > > i = I, 2,3 denote the rows of 

\ curl 4___ / 


On the one hand, considering the variation 5e^ € 0^(0, T) of the left hand side of 



Existence results in rate-independent gradient plasticity 


13 


Additive split of distortion: 
Equilibrium: 

Free energy: 

Yield condition: 
where 


Vu = e + p, e®:=syme, := symp 
Div cj + / = 0 with a = C.e® 

jCurlePp + ^/iA:2 |7p 
(j){T,E,g) ■■= |devsymS£;| -\- g - ay < 0 
EE:=a + = -pL^ Curl Curl 


g = -hk2i 


Dissipation inequality: 

Dissipation function: 

Flow law in primal form: 

flow law in dual form: 

KKT conditions: 

Boundary conditions fore^: 
Function space for : 


[ [{^E,£^) +9i] dx>0 
JQ 

ay\q\ 


^(9 ,6 := 


oo otherwise 
{EE,g) € dV{eP,j) 

devsymSE ^ = A = 1^ 

I dev sym L e \ 

A > 0, (j){EE,g) <0, X(j){EE,g) = 0 

X n = 0 on r, (Curl e^) x n = 0 on dD \ T 
e^(t, •) G H(Curl; n, Sym(3)) 


Table 5; The irrotational version of [17] with isotropic hardening. The plastic distortion 
itself does not appear, only = sym p remains in the model. 


(3.1) with respect to the plastic strain variable we get 

[ p L^l Curl {eP + t6£P)\‘^ dx\ = pL^ [ (Curl Curl (5e^) dx 


dt 2 


t=o 


p y' / (curl £^, curl fef) dx = 4 pL^ V' / (axl skew Vef, axl skew Vfef) dx 
i=i i=i 

3 „ 3 „ 

/ (skew Vef, skew Vfef)dx = J^ / {2pLl skew , Vdef) dx 


2=1 


2=1 


m^, 2nd order tensor 
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= “IZ / (Divmf,fef)dx + ^ f (mf.n,fef)(ia 
i=i i=i J9n 

3 » 3 „ 

= — ^ / (Divm^, 6e^) dx + ^ / ((axlm^) x n, da 

i=i i=i 

3 „ 3 „ 

= — ^ / (Divmf,fef) dx — ^ / (axlmf, (&f) x n) da 

i=l i=l 

= — f {Div xxF, 5e^) dx — j (axlm?, (de^) x n) da, 

dn dan 


where the third order tensor € j^sxsxs -g (Jef^j^ed by 


mP = (m?,m^,mf)'^. 


(3.2) 


(3.3) 


On the other hand, we obtain 


f Curl (e^ + tde^)pdx [ (Curl e^. Curl de^) dx 

dt 2 dn i=o dn 

= f (Curl Curl s^, de^) dx + fiL‘1 ^ f (def x [Curl £^]i, n) da (3.4) 
Jn dan 

= IJ,L^ f (Curl Curl s^, de^) dx — fiL^ 'S^ f ([Curl e^]i, x n) da. 
dn dan 

In view of (3.2) and (3.4), we obtain 

laL^ f (Curl Curl de^) dx = — f (Divm^, de^) dx, (3.5) 

J Q. J Q. 

for all 6eP G Ho(Curl; fl; P), i.e., for def x n = 0. 

Since we can assume that is trace free symmetric, so is 5eP and we may equiva¬ 
lently write 

fiL^ f (dev sym Curl Curl e^, de^) dx = — f (Divm^, de^) dx . 
dn dn 

Thus, we get that 

—Div = fiL^ dev sym Curl Curl . 

6K3X3X3 ' 37^ ' 


tP := dev a + Div = dev a — dev sym Curl Curl 


= dev sym E e 


Set 


(3.6) 
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Hence, = devsymS^;. 

Notice that the second order tensor Divm^ is trace free. In fact, from the bracket 
(m^, Ve^), it holds that (we may assume) (since is symmetric) and we 

nay also assume that 

<fc = 0. (3.7) 

Now, it is clear that tr(Divm^) = muk^k = 0 from (3.7) and hence Divm^ is symmetric 
and trace free. 

Find a summary of this model in Table 6. 

Let us now repeat the formulation of the model with spin in more details, in its 
dual and primal setting for the paper to be rather self-contained. 


4 The model with linear kinematical hardening and plas¬ 
tic spin 

4.1 Strong formulation 

The balance equation. The conventional macroscopic force balance leads to the 
equation of equilibrium 

div (7 + / = 0 . (4.1) 

Constitutive equations. The constitutive equations are obtained from a free energy 
imbalance together with a flow law that characterizes plastic behaviour. Since the 
model under study involves plastic spin, we consider an additive decomposition of the 
displacement gradient Vu into elastic and plastic components e and p as mentioned in 
the notational section, so that 

(4.2) 

Ti-i(Curlp) (4.3) 

elastic energy defect energy (GND) 

+ gi.(£ ) 

linear kinematical hardening energy 

4'e'‘(e) := ^(e^C.e^), T(,“i(Curl p) := \ Curl and 

^kiniP) ■= Idevsympp . 


Vu = e + p 

We consider here a free energy of the form 

T(Vii,p, Curl p) ; = 'I'e°(e) + 


where 
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Additive split of strain: 
Equilibrium: 

Microforce balance: 
where 

sym Vu = £ = + £P , G Sym (3) 

Div cj + / = 0 with a = C.e® 
dev a = t'p — Div 

T^: microstress ( 2 “^ order) 
m^: micropolar stress (3'^'^ order) 
tr(DivmP) = 0 

t'p = devfj + Divm^ G s[(3) n Sym (3) 

Free energy: 

+ \p.L‘l jCurlePp + i/iA: 2 | 7 p 

Yield condition: 

■= \t^\ +g-cfy<b 

where 

mP = mfnerg 


pLl (Curl eP, CurleP) = (mP, Ve^) 

Dissipation inequality: 

[ [(rP,ffP) + gf] dx>0, g = -pk 2 y 

JQ 

Dissipation function: 

1 kl if \q\ < C, 

' ’ ^ ( +00 otherwise 

Flow law in primal form: 

{TP,g) G dV{£P,j) 

Flow law in dual form: 

tP tP 

= A-= ^, 7 = A = |e^’| 

ay-g JrP] 

KKT conditions: 

A > 0, (j){TP,g)<0, X(j){TP,g) = 0 

Boundary conditions fors^: 
Function space for : 

Length scale: 

eP X n = 0 on dQ 

£P{t,-) G H(Curl; D, Sym (3)) 

energetic Lc 


Table 6 ; The irrotational model by Gurtin and Anand [21] with no dissipative length 
scale i.e. = 0. Since has been identified with devsymS^; in (3.6), the model 

coincides with the irrotational version of [17]. 


Lc is the energetic length scale and ki is the dimensionless hardening modulus. 
The defect energy is related to geometrically necessary dislocations (GNDs) and the 
Burger’s vector. 

The local free-energy imbalance states that 


T - {a,e) - {a,p) < 0 . 


(4.4) 
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Now we expand the first term, substitute (4.3) and get 

(Ce® — (T, e®) — ((7,p} + pLl (Curlp, Curlp) + pki (devsymp,p) < 0, (4.5) 

which, using arguments from thermodynamics gives the elasticity relation 

a = = 2p sym(Vu — p) + A tr(V — p)l (4.6) 


and the reduced dissipation inequality 

— {a,p) -\- phi (Curl p. Curl p) + pki (devsymp,p) < 0. (4-7) 

Now we integrate (4.7) over 14 and get 


0 > 


/ + p (Curl p. Curl p) pki (dev symp,p) 

Jn L 

[ [“(<7)P) + pL^ (Curl Curlp,p) -\- pki (devsymp,p) 


dx 


d ,• 


+ ^div(pL2 _pi X (Curlp)* 


2=1 


dx. 


(4.8) 


Using the divergence theorem we obtain 


In 


[(—cj + pL^ Curl Curl p,p) + p/ci (devsymp,p)] dx 
p lI (p* X (Curl p)*, n)da <0 . 


2=1 




(4.9) 


In order to obtain a dissipation inequality in the spirit of classical plasticity, we 
assume that the infinitesimal plastic distortion p satisfies the so-called linearized insu¬ 
lation condition 

3 ^ , 

pLl{— p^ X {Curl py,n)da = 0. (4.10) 

^ Jan dt 

This condition is satisfied if we assume for instance that the boundary is a perfect 
conductor. This means that the tangential component of p vanishes on did. In the 
context of dislocation dynamics these conditions express the requirement that there is 
no flux of the Burgers vector across a hard boundary. Curtin and Anand [21] introduce 
the following different types of boundary conditions for the plastic distortion 


dtp X n|riiard = 0 ’’micro-hard” (perfect conductor) 
^iPlrhard = 0 ” hard-slip” 

Curl p X n|rhard ~ 0 ’’micro-free” . 


(4.11) 


We specify a sufficient condition for the micro-hard boundary condition, namely 


P X ^Irhard = 0 


(4.12) 
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and assume for simplicity only Phard = d^l = P. Note that this boundary condition 
constrains the plastic slip in tangential direction only, which is what we expect to 
happen at Phard- 

Under (4.10), we then obtain the dissipation inequality 

[ {a + E^^,, + El^,p)dx>0, (4.13) 

Jn 

where 

^curi Curl Curl p and := —devsymp. 


The flow law. We consider a yield function defined for every Se := 
by 

00(Se) := I dev Se| - o-y 


cj + S 


lin 

curl 


+ s 


lin 

kin 


(4.14) 


Here Uy is the yield stress of the material. So the set of admissible (elastic) generalized 
stresses is 

/Co := {Se = a + I 0 o(Se) < o} . (4.15) 

The maximum dissipation principle gives the normality law 


P€N,Co{^e)) (4.16) 

where N!Cq{Ee) denotes the normal cone to /Co at Se, which is the set of generalised 
strain rates p that satisfy 


(S — Ee,p) < 0 for all S G /Co . 


(4.17) 


Notice that Nic^ = dXo where Xq denotes the indicator function of the set /Co and dXo 
denotes the subdifferential of the function Xq- 
Whenever the yield surface dJCo is smooth at then 


P G Nico{Ee)) 


3 A such that p 


dev Se 
dev Se| 


with the Karush-Kuhn Tucker conditions; A > 0, 0 (Se) < 0 and A0 (Se) = 0. 
Using convex analysis (Legendre-transformation) we find that 


p G 0Xo(Se) ^ Se G dXUp), 


(4.18) 


where Xq is the Fenchel-Legendre dual of the function Xo denoted in this context by 
"D, the one-homogeneous dissipation function for rate-independent processes. That is, 

V{q) = sup {(a + 1 + , q) \ (^o{a + ) < o} = Uy |g| . (4.19) 

We get from the definition of the subdifferential (Se G c/Xo( 0)) that. 


Pio) P P{p) + (51e, Q — p) for any q. 


(4.20) 
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That is, 

Vig) > V{p) + {a + 1 + , q - i>) for any q. (4.21) 


Strong formulation of the model. To summarize, we have obtained the following 
strong formulation for the model of infinitesimal gradient plasticity with kinematic 
hardening and plastic spin. The goal is to find: 

(i) the displacement u € H^(0, T; HQ(n, T,R^)), 

(ii) the infinitesimal plastic distortion p with symp G H^(0, T; L^(n,s[(3))), Curl p G 
Hi(0,r;L2(fl,R3x3)) and Curl Curl p G 



Table 7: The model with linear kinematical hardening and plastic spin. The boundary 
condition on p necessitates at least p G H(Curl; fl, R^^^), However, whether this is the 
case will only be proven at the end of the paper. 
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4.2 Weak formulation of the model 

Assume that the strong formulation has a solution Let v G with 

U|r = 0. Multiply the equilibrium equation with v — u and integrate in space to get 




Using the symmetry of the stress tensor a and the elasticity relation we get 


(4.22) 


(C. sym(Vu — p), sym(Vu — Vu))dx 


/ f{v — u)dx . 

Jn 


(4.23) 


Now, we take any q G C'°°(n,sl(3)) such that g x n = 0 on F and we integrate (4.21) 
over n, integrate by parts the term with Curl Curl using the boundary conditions 


(g — p) X n = 0 on r 


and get 

[ V^^^q)dx > 
Jn 

> 


> 


VQ^'^{p)dx+ / (cj + S“,i + S])“,g-p)dx 

I Jn 

Po^^ip) dx + / (C. sym(Vu — p), sym(g — p)) dx 
I Jn 

[ {pLl Curl Curl p + pki devsymp, q — p) dx 


Po{p)dx+ / (C. sym(Vu — p), sym(g — p)) dx 


(4.24) 


— pL^ / (Curl p, Curl (g — p)) dx — p fci / (symp, g — p) dx . 
Jn Jn 

Adding (4.24) to the weak formulation of the equilibrium in (4.23), we get that 


In L 


(C. sym(Vu — p), sym(Vu — q) — sym(Vh — p)) + /iL^(Curl p. Curl {q — p)) 


+pA:i(symp, symg — symp) dx + / I?g^“(g)dx— / Pg™(p)dx 


-.kin/ 


> / f{v — u)dx y{v,q) 

Jn 


(4.25) 


4.3 Existence result for the new formulation 

To prove the existence result for the weak formulation (4.25), we follow the abstract 
machinery developed by Han and Reddy in [24] for mathematical problems in classical 
plasticity and used for instance in Djoko et al. [14], Reddy et al. [46], Neff et al. [34], 
Ebobisse-Neff [17] for models of gradient plasticity. To this aim, (4.25) is written as 
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the variational inequality of the second kind: find w = {u,p) € H^(0,T;Z) such that 
t(;(0) = 0 and 

a{w, z — w) + jo{z) — jo{w) > {i,z — w) for every z £ Z and for a.e. t G [0, T] , (4.26) 
where Z is a suitable Hilbert space to be constructed later, 


a(w,z) := / (C.(sym(Vu — p)), sym(Vu — g)) +//L^(Curl p, Curl g) (4.27) 

Jfi L 

+p ki (symp, sym g) 1 dx , 


Mz) := [ 75o“(g) dx , (4.28) 

Jn 

{e,z) := [ fvdx, (4.29) 

Jn 

for w = {u,p) and z = {v,q) in Z. 

The Hilbert space Z is constructed in such a way that the functionals a, jo and I satisfy 
the assumptions in the abstract result in [24, Theorem 7.3]. The key issue here is the 
coercivity of the bilinear form a on Z. From the structure of the bilinear form a and 
the functional jo, a natural attempt for the space of infinitesimal plastic distortions, is 
to consider the closure Mgym (Curl, H,r;s[(3)) of the linear subspace 

{q € I trg = 0, g x n = 0 on T} 

with respect to the norm 


lkllsym,curi — l|sym911^2 + ||Curlg ||^2 . (4.30) 

Motivated by the well-posedness question for our model [34, 17], Neff et al. [37, 38, 
39, 40], derived a new inequality extending Korn’s inequality to incompatible tensor 
fields, namely there exist a constant C'(H) > 0 such that 

Vp€ H(Curl; 12, p x n|r = 0 : (4.31) 

\\p\\l^{ ni < C{n) ( II symp||j,2(^)^ + || Curl ph^^n)^ ) • 

plastic distortion plastic strain dislocation density 


Here, T C 912 with full two-dimensional surface measure and the domain 12 needs to 
be sliceable, i.e. cuttable into finitely many simply connected subdomains with Lip- 
schitz boundaries. The inequality (4.31) expresses the important fact that controlling 
the plastic strain symp and the dislocation density Curl p in L^(12) gives a control 
of the full plastic distortion p in L^(12) provided the correct boundary conditions are 
specihed: namely the micro-hard boundary condition. Since in the sequel we assume 
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that tr{p) = 0 (plastic incompressibility) the quadratic terms in the thermodynamic 
potential provide a control of the right hand side in (4.31). So, setting: 

V = Hj(54,r,M3) = {u € H^(f2,R3) |u|r = 0}, (4.32) 

Q = Ho(Curl; 12, r,sl(3)) = Hsym (Curl; f2,r;sl(3)), (4.33) 

Z = V X Q, (4.34) 

equipped with the norms 

:= l|Vu ||^2 , II^IIq := ||symg ||^2 + ||Curlg ||^2 , (4.35) 

ll^llz := ll^lly + Ikllg for z = (u, g) G Z . (4.36) 

Let us show that the bilinear form a is coercive on Z. Let therefore z = (v, q) G Z. 
a{z,z) > mo ||sym Vu — symg|| 2 +/iLc ll^url g ||2 +/Ufci ||symg ||2 (from (2.2)) 

= mo ||sym Vu||2 + ||symg||2 — 2(sym Vu,symp) 

+ pLl IlCurl g||2+ ^A:i ||symg||2 


> mo 


||sym Vu ||2 + ||symg ||2 - 6 ||sym Vu ||2 - ^ ||symg ||2 

||Curl q\\2 + pki ||symg ||2 (using Young’s inequality) 


= mo(l — 6*) ||sym Vu ||2 + 

mo 


mo(l + pki 


\\sy^Q\\l+ fLI IlCurl q\\l. 


So, choosing 6 such that 


mo + pki 

find a positive constant C{mo, p, ki,Lc, 12) > 0 such that 


< 6 < 1 and using Korn’s first inequality, we 


a{z,z) > C \\v\\y + ||symg ||2 + ||Curl q\\l 


= C'lkl 


Vz = (u, g) G Z , 


which proves the coercivity of our bilinear form and the inequality (4.31) shows the 
equivalence Q = Hsyjn(Curl,12,r;sl(3)). 


5 The Gurtin-Anand model with linear kinematical hard¬ 
ening: purely energetic version 


Constitutive equations. 

Vu = e + p symVu = e® + e^. (5.1) 

We consider here a free energy of the form 

^(e^eP,Curle^’) ; = ^e“(e") + ^|i“i(Curl e^) (5.2) 

elastic energy defect energy (GND) 



linear kinematical hardening energy 
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where 

^!!°i(CurUP):=iMic|Curle?’|2 and 
^kTn(e^) := ^pki\deYeP\^ . 

Following the development in section 4, the free energy imbalance taking into account 
the boundary condition of the plastic strain variable 

X n|r = 0 


leads to the dissipation inequality 


{T,E,e^) dx > 0 , 


(5.3) 


where 


He ■= cr — pki dev Curl Curl . 

The flow law. The set of generalized stresses is 

JC := {E € Sym (3) | (j){Tj) := \ dev S| — Uy < 0} . 
Hence, following [17], we get the flow in dual form 

eP G N/cCEe) 


(5.4) 


where N/cCEe) denotes the normal cone to fC at Ee, which in case of smoothness reads 
as 

dev Ee 

with A > 0, (()(S_b) < 0, X(P(Ee) = 0. 

The flow law in its primal formulation reads as 

Se G dV{eP) . 

That is, 

^{q) > V{eP){'EE,q - eP) VgGSym(3), (5.5) 

where "D is the dissipation function defined as 

^( 9 ) := {(5^> 9 ) I S G /C} = fjy \q\ Vg G Sym(3). 

Weak formulation of the model. Now arguing as in Section 4 and also as in 
Ebobisse-Neff [17, Section 3], we obtain a weak formulation of the model in the form 
of a variational inequality 

(C.(sym Vrt — e^), (sym Vu — q) — (sym Vh — e^)) + pL^{Curl e^. Curl {q — e^)) 
+/i/ci(e^, g — e^) dx-\- f V{q)dx— f 'D{e^)dx 

^ Ju Jo 

> [ f{v — u)dx V(u,g). (5.6) 
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The existence and uniqueness result for the variational inequality is easily obtained 
in the spaces 

u G H^(o,r; r, R3)), 

sP € H^(0,r; Ho(Curl,r;Sym(3) nsl(3))), 
as in Section 4 through [24, Theorem 7.3], following the coercivity on the space 
Z := r, X Ho(Curl, T; Sym(3) n s[(3)), 


of the bilinear form 


sl{w,z):= / (C.(sym Vu — p), (sym Vu — g)) + ;U L^(Curl p. Curl g) + /ci(p, (7) 


dx 


for every w = {u,p), z = {v,q) in Z. 

Note that since is already trace-free and symmetric, the coercivity of the bilinear 
form a does not the new Korn’s type inequality in [37, 38, 39, 40], and in (4.31). 


6 The infinitesimal elastic micromorphic model 


The same total energy 

£{u,p) = , 


in 


(C. sym(Vu — p), sym(Vtt — p)) 


pki 


pLl 


+^-^| devsympj^ H-Curl p|^ — (/, u) 


dx 


( 6 . 1 ) 


2 ' 2 
is the starting point for a two-field minimization formulation 

S{u,p) min. w.r.t {u,p), 

in the sense of a micromorphic model ([35, 36]). 

The relation of (6.1) to our plasticity formulation (1.2) is that in (6.1) the micromorphic 
distortion p is determined directly by a global energy minimization instead of a plastic 
flow rule. The microbalance equation is obtained as follows. The first variation of (6.1) 
with respect to p gives 

/ [(C. sym(Vu — p), sym (5p)-|-p fci(devsym p, 5p)-|-/i (Curl Curl p,Sp)]dx 

Jn 

= / [(sym C. sym(Vu — p), (5p)-|-(p/ci, devsym p, (5p) 

Jn 

+ (p L‘1 Curl Curl p, 5p)] dx = 0 . 

The ’’microbalance” is then of the form 

Cauchy stress a 
' -' 

p Curl Curl p = sym C. sym(Vu — p) —p ki dev sym p. 
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The well-posedness of such a model has been shown in Neff et al. in [36]. Hence, we 
get in this model 

p Curl Curl p = a — pki dev sym p 


or 


t) = a — pki dev sym p — pL^ Curl Curl p 




with 


p X n|r = 0 and (Curl p) x = 0 ) 

instead of a dual flow law p G BX^Ee) in plasticity. 


7 Conclusion 

The development of the model with plastic spin is straightforward and involves only 
the addition of a quadratic defect energy. The boundary conditions on the plastic dis¬ 
tortion are consistent both from the physical and the mathematical point of view. The 
departure from classical plasticity is minimal. Choosing a symmetric local kinematical 
backstress evolution necessitates to use a new Korn’s type inequality for incompatible 
plastic distortions. Contrary to the presented alternative models, in which the energetic 
length scale L^. has only a ’’passive” role in that necessary estimates are already ob¬ 
tained from the dissipative length scale i, in this model it is only the interplay between 
the energetic length scale and the symmetric local backstress which makes the problem 
well-posed. By identifying the irrotational Gurtin-Anand model with only energetic 
length scale as a special limit of our model with spin, we have been able to provide an 
existence theorem for that model for both the isotropic hardening case ([17]), as well as 
the local backstress case (this paper with the same considerations as in [17]). Moreover, 
our derivation of the model avoids the introduction of certain additional ’’micro force 
balances”. Let us also mention that, the introduction of the irrotationality constraint 
appears, in our general framework with spin to be neither advantageous nor necessary, 
but simplifies the analysis considerably. 

It remains to be seen if, in the dual formulation of our model with spin one may 
consider isotropic hardening driven by a symmetrized measure of accumulated plastic 
straining 

7 = |symp|. (7.1) 

This would be conceptionally pleasing since kinematical hardening could then exclu¬ 
sively be related to the GND-distribution via the energetic length scale Lc (and as¬ 
suming ki = 0) while the SSD-distribution would be described by the accumulated 
plastic straining. We need to remark that (7.1) does not seem to satisfy the additional 
assumption of maximal dissipation, making it unsuitable to be considered in the primal 
formulation. However, it is well-established that the equivalence of the primal and dual 
formulation is not satisfied in general for gradient plasticity. 
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